NARROW ESCAPE, part II: The circular disk 



A. Singer * Z. SchuW, D. Holcman^ 
February 7, 2008 



Abstract 

We consider Brownian motion in a circular disk f2, whose boundary 
<9f2 is reflecting, except for a small arc, d£l a , which is absorbing. As 
e = |9f] a |/|<9ri| decreases to zero the mean time to absorption in dfl a , 
denoted Et, becomes infinite. The narrow escape problem is to find an 
asymptotic expansion of Et for £< 1. We find the first two terms in 
the expansion and an estimate of the error. The results are extended in 
a straightforward manner to planar domains and two-dimensional Rie- 
mannian manifolds that can be mapped conformally onto the disk. Our 
results improve the previously derived expansion for a general smooth 



domain, Et = 

Dir 



io g i + o(i; 



(D is the diffusion coefficient) in the 



case of a circular disk. We find that the mean first passage time from the 



R 2 

center of the disk is E[t \ x(0) = 0] = -jj 



log^ + log2 + i + 0(e) 



The second term in the expansion is needed in real life applications, 

such as trafficking of receptors on neuronal spines, because log — is not 

e 

necessarily large, even when e is small. We also find the singular be- 
havior of the probability flux profile into <9f2 a at the endpoints of dO, a , 
and find the value of the flux near the center of the window. 



1 Introduction 

The expected lifetime of a Brownian motion in a bounded domain, whose 
boundary is reflecting, except for a small absorbing portion, increases indefi- 
nitely as the absorbing part shrinks to zero. The narrow escape problem is to 
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find an asymptotic expansion of the expected lifetime of the Brownian motion 
in this limit. The narrow escape problem in three dimensions has been studied 
in the first paper of this series £Q , where is was converted to a mixed Dirichlet- 
Neumann boundary value problem for the Poisson equation in the domain. 
This is a well known problem of classical electrostatics (e.g., the electrified 
disk problem j2]), elasticity (punch problems), diffusion and conductance the- 
ory, hydrodynamics, and acoustics 0-0- It dates back to Helmholtz [S] and 
Lord Rayleigh P and has been extensively studied in the literature for special 
geometries. 

The study of the two-dimensional narrow escape problem began in 
in the context of receptor trafficking on biological membranes ^T], where a 
leading order expansion of the expected lifetime was constructed for a general 
smooth planar domain. In this paper we present a thorough analysis of the 
narrow escape problem for the circular disk and note that our calculations 
apply in a straightforward manner to any simply connected domain in the 
plane that can be mapped conformally onto the disk. According to Riemann's 
mapping theorem ^2] 5 this covers all simply connected planar domains whose 
boundary contains at least one point. The same conclusion holds for the 
narrow escape problem on two-dimensional Riemannian manifolds that are 
conformally equivalent to a circular disk. The biological problem of receptor 
trafficking on membranes is locally planar, but globally it is a problem on a 
Riemannian manifold. The narrow escape problem of non-smooth domains 
that contain corners or cusp points at their boundary is treated in the third 
part of this series > where the conformal mapping method is demonstrated. 

The specific mathematical problem can be formulated as follows. A Brown- 
ian particle diffuses freely in a disk Q, whose boundary dQ is reflecting, except 
for a small absorbing arc dQ a . The ratio between the arclength of the absorb- 
ing boundary and the arclength of the entire boundary is a small parameter 



\dn a \ 

\dtt\ 



<C 1. 



The mean first passage time to dQ a , denoted Et, becomes infinite as e — > 0. 
The asymptotic expansion of Et for e < 1 was considered for the particular 
case when dfl a is a disjoint component of dfl in j!4[ and references therein]. 
This case differs from the case at hand in that the absorption probability flux 
density in the former is regular, while in the latter it is singular. It was shown 
in [TU] that Et for the narrow escape problem in a general planar domain Q 
has the asymptotic form 



Et 



\n\ 



io g i + o(i; 



'i.r, 



where \Q\ is the area of Q, and D is the diffusion coefficient. This leading order 
asymptotics has the drawback that loge can be 0(1) when e < 1. Thus the 



2 



second term in the expansion is needed. For the particular case of a circular 
disk an approximate value for the correction was given in ^Uj- in contrast, the 
asymptotics of Et for a three dimensional ball of radius R with an absorbing 
window of radius sR is 



Et 



\n\ 



ADeR 



[l + 0(eloge)] 



so the leading order term is much larger than the correction term if e is small. 
The difference in the asymptotic form of Et stems from the different singular- 
ities of the Neumann function in two and three dimensions: it is logarithmic 
in two dimensions and has a pole in three dimensions. 

Our computations are based on the mixed boundary value techniques of [3] . 
They reveal the singularity of the absorption flux in the absorbing arc dQ a . 



Specifically, the singularity is (e 2 



"V 2 , where s is the (dimensionless) 



arclength measured from the center of d£l a , and attains the values s = ±e at 
the endpoints. 

The exit time vanishes at the absorbing boundary, and is small near the 
absorbing boundary, but it attains large and almost constant values of order 

log - inside the domain. We show that this "jump" occurs in a small bound- 
er 

ary layer of size O log —J . We calculate the average exit time, where the 

averaging is against a uniform initial distribution in the disk, the time to exit 
from the center, and the maximum mean exit time, attained at the antipodal 
point to the center of the absorbing window. 

The mean first passage time (MFPT) from the center of the disk is 



E[t\x(0) =0] 



D 



log^ + log2 + ^ + 0(e) 



;i-2) 



the MFPT, averaged with respect to an initial uniform distribution in the disk 
is 

1 



Et 



D 



log - + log 2 

e 



0(e) 



;i.3) 



and the maximal value of the MFPT is attained on the circumference, at the 
antipodal point to the center of the hole, 

= 0] 



max Mr I x] = E\r I r = 1, 
xen 



IP 
D 



log - + 21og2 + 0(e) 



;i-4) 



The boundary layer analysis of Et can be applied to the ap- 
proximation of the first eigenfunction and eigenvalue of the mixed 
Neumann-Dirichlet boundary value problem with a small Dirichlet 
window on the boundary. This problem arises in the construction of 
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the first eigenfunction and eigenvalue of the Neumann problem in a 
domain that consists of two domains (e.g., circular disks) connected 
by a narrow channel [15j, [16J. 

Specifically, it is easy to see that 



Et 



^ 1 1 

— ' Ar, An 



1.5) 



ra=0 



where < Ao <C Ai < • ■ ■ are the eigenvalues of the mixed problem 
and the MFPT is also averaged with respect to the initial point. 
The first eigenfunction w n of the mixed problem is differs from the 
first eigenfunction of the Neumann problem, which is Vq — 1, only 
in a boundary layer about the small window. Thus u is a small 
perturbation (in L 2 norm) of Vq = 1. It follow that w /An differs from 
Et only in the boundary layer. 

2 Solution of a mixed boundary value problem 

In non-dimensional variables the narrow escape problem concerns Brownian 
motion inside the unit disk, whose boundary is reflecting but for a small ab- 
sorbing arc of length 2e (see FigJTJ). In polar coordinates x = (r, 9) the MFPT 

v(r,9) = E[t\x(0) = (r,9)], 

is the solution to the mixed Neumann-Dirichlet inhomogeneous boundary value 
problem (see, e.g. [T7j ) 

Av(r,6) = -1, r<l, for < 9 < 2vr, 



v(r, 6) 
dv{r, 9) 



Or 



r=l 



r=l 



0, for |^ — 7r| < e, 



0, for \6-ir\)e, 



which is reduced by the substitution 



u = v 



1-r 2 
4 



(2.1) 



(2.2) 



to the mixed Neumann-Dirichlet problem for the Laplace equation 



Au(r, 



0, for r < 1, < 9 < 2vr, 



u(r, 9) 
du(r, 9) 



Or 



r=l 



r=l 



0, for \9 — 7r| < £, 



-, for \e-ir\)e. 



(2.3) 
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We adapt the method of [3] to the solution of (|2.3|) . Separation of variables 
suggests that 

CO 

Oo 

2 



u(r, 9) = y + a nr n cos n6, (2.4) 

n=l 

where the coefficients {a n } are to be determined by the boundary conditions 

oo 

u(r,9) = — + ^ a n cosn6 = 0, for ir — e < 9 < n, (2.5) 

r=l 2 n=l 

- ^ na n cos n0 = |, for < 9 < tt - e. (2.6) 



dr 



r=l 



n=l 



We identify this problem with problem (5.4.4) in where general functions 
appear on the right hand sides of equations ()2.5|) (j2.6|) . Due to the invertibility 
of Abel's integral operator, the equation 



+ a n cos n9 = cos -6* 



n=l 



^l(i) d£ 



a/ cos 9 — cos £ 



(2.7) 



defines hi(t) uniquely for < £ < n — e. The coefficients are given by 



cos n9 cos -9 d6 
2 



hx{t) dt 



h x {t) dt 



y/ cos 9 — cos £ 

cos (n + ~) 9 + cos (n — |) # 
V cos 6* — cos £ 



d0. 



The integral 



P„(cosu) = 



y/2 f u cos(n + ±)9 



de. 



X8) 



;2.9) 



7T Jq V COS 9 — COS M 

is Mehler's integral representation of representation of the Legendre polynomial 
[TH] . It follows that 



* v/2 Jo 

for n > 0, and 

£ii(£)rf£ 
Integration of ()2.6|) gives 



/ii(£)[P n (cos£) + P n _i(cos£)] dt, 



(2.10) 



2 

a = - 

7T 



cos -9 

2 d9 = V2 I h x {t)dt. (2.11) 



o Vcos# — cos£ 



oo . 

a ra sin n9 = -9, for < 9 < n — e. 



;2.i2) 



n=l 



Changing the order of summation and integration yields 

/ hx{t)—= V[P n (cost) + P n _i (cost)] sin n#dt = -9. (2.13) 
^ v2 n=1 2 

Using P eq.(2.6.31)], 

1 ^ r , , cos yH(0-t) 

-= V[P n (cost) + P n _i(cost)] sinnfl = - 2 (2.14) 
V 2 , V cos £ — cos 

v n=l v 

we obtain 

\ - 1W 1-, for 0<^<tt-£. (2.15) 

Jo Vcost-cos^ 2 cos 2^ 

The solution of the Abel-type integral equation ()2.15|) is given by 

u 

h x {t) = -- 2 du. (2.16) 

71 di J y/ cos u — cos t 

Together with ()2.11j) this gives 

u 

ny rir-e u sin — 
a = — 2 = du. (2.17) 

71 Jo v cos M + cos e 

We expect the function w(r, #), closely related to the MFPT, to be almost 
constant in the disk, except for a boundary layer near the absorbing arc. The 
value of this constant is a , because all other terms of expansion (|2.4j) are 
oscillatory. 

2.1 Small e asymptotics 

The results of the previous section are independent of the value of e. Here we 
find the asymptotic of a for e <C 1. Substituting 



/cosu + coee 
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in the integral (|2.17|) yields 



4 , cos ( e /2) arccos J s 2 + sin 2 - 
a = — / = ds 

71 -'o / 2 , • 2 £ 

■ 1 s z + sin - 



-cos(e/2) 1 a /-cos(e/2) arcsm A / S + Sm 9 

2/ da - - / , V cis 



>2 I mr,2 5 



2 



£ 4 / .cos( e /2) arcsm Ws 2 + sin 2 | 
2 log ( 1 + cos - ) — 2 log sin / = — (is 

V 27 2 7T J Q / . 2 £ 

V s +Sm 2 

= -21og^ + 21og2-l f'^^ds + Oie) 

2 7T J S 

= -2 log I + 0(e), (2.19) 
Z" 1 arcsin s -k 

because / ds = —log 2. The substitution (12.181) turns out to be 

Jo s 2 
extremely useful in evaluating the integrals appearing here. 

2.2 Expected lifetime 

Now, that we have the asymptotic expansion of a (eq. ()2.19|) ). the evaluation 
of expected lifetime (MFPT to the absorbing boundary dQ a ) becomes possible. 
Setting r = in equations ()2.2|) and (|2.4|1 . we obtain the expression (|1.2|) for 
MFPT from the center of the disk. 

Averaging (|2.H|) with respect to a uniform initial distribution in Q gives 



j /-27T rl 



71 



-Er — — I dO I v(r,9)rdr 









r 3 


f 


(a -\ 




r 


h 


~2y 


2 _ 



dr 



as asserted in eq. (jl.3|) . 

The maximal value of the MFPT is attained at the point r = 1, 9 = 0, 

which is antipodal to the center of the absorbing arc. At this point — = 0, 

__ ad 
as can be seen by differentiating expansion (J2.4|) term by term. Setting r = 1 
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and 6 = 0, we find that 

oo 

Vmax = u(l, 0) = y + a n . (2.21) 

n=l 

The evaluation of the maximal exit time is not as straightforward as the 
previous evaluated MFPTs, because one needs to calculate the infinite sum in 
f!2.21|) . This calculation is done in Appendix where we find feq. (jA.12|) ) 

Vmax = log ^ + 2 log 2 + 0(e), 

as asserted in equation (jl.4|) . 



2.3 Boundary layers 

We see that the maximal exit time is only v max — v center = log 2 — — = 

.4431471806 . . . longer than its value at the center of the disk. In other words, 
the variance along the radius = 0, < r < 1 is very small. However, in 
the opposite direction 6 = it, < r < 1, we expect a much different behavior. 

In particular, the MFPT is decreasing from a value of v cen t er ~ log - at the 

e 

center of the disk to v(l, tt) = at the center of d£l a . The calculation of the 
exit time 

v ray {r) = v(r, 9 = n) = + ^ + ]T a n (-r) n , (2.22) 



n=l 



is similar to that of the maximal exit time and is done in Appendix El For 
e 1 and 1 — r ^> y/e, we find the asymptotic form (eq. ([B.8[l ) 



v ray (r) = - log | + 2 log(l - r) + - log(l + r 2 ) + q(r) +0(e), (2.23) 

where q(r) is a smooth function in the interval [0, 1] (eqs. ([B.6j) - ([B.7)l ). Clearly, 
this asymptotic expansion does not hold all the way through to the absorbing 
arc at r = 1, where the boundary condition requires v ray (r = 1) =0. Instead, 

the boundary condition is almost satisfied at r = 1 



VrayU-Jl ) = -log ~ + 2\og Ufj + 0(s) = 0(e), (2.24) 

In other words, the asymptotic series (|2.23|) is the outer expansion [T5] . 
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We proceed to construct the boundary layer for 1 — r y/e. Setting 
5 = 1 — r, we have the identities 



1-r 2 



2 4 ' 



l-2rcos£ + r 2 = 4 sin 2 - (1 - 5) + 5 2 . 



The exact form of the MFPT along the ray, eq. (|B.3J) . gives the expansion 

Vray(S) = £ + ( 2 -2f> ) 



9 A • fc 

z 4 sin - 
2 



5 



7r sm - J o 
2 



r „ M;/2) arccos J s 2 + sin 2 - s 2 ds 



+ sm - I 

2 / 



+ 



<5 5 



Evaluating the integral in eq. (|2.25|) . 



<nM 2) arccos ^/s 2 + sin 2 | s 2 



s 2 + sin 2 



2, 



7T 
2 



log sin ^ + cos ^ — log ^1 + cos ^ ) log 2 



we obtain the boundary layer structure 



Vray(S) 



5 7 



- + 8, - . 



+ 0(e), (2.26) 



[2.27) 



In particular, setting 5q = — £ l°g 2 yields 



*Va#o) = - log- + 0(elog 2 e), 



(2.28) 



which is the value of the outer solution. We conclude that the width of the 
boundary layer is O I e log - ) . Furthermore, the flux at the center of the hole 



is given by 



A ux center 



dv r 
dr 



ay 



du 



ray 



r=l 



85 



0(1). 



(2.29) 



5=0 
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2.4 Flux profile 

Next, we calculate the profile of the flux on the absorbing arc. Differentiating 
expansion (|2.4|) gives the flux as 



du(r, 9) 



r=i dr 



1 1 °° 

- = h na„ cos n9, (2.30) 



r=1 2 » - 1 



9r 

for 7r — e < 9 < 7r. Using equation (J2.10)) for the coefficients, we have 

/(0) = — + ~r fatf) dtJ2n[P n (cost) +P n -i (cost)] cos n6 

1 Id r 7r ^ £ 
= ~2 + ~^d9j M A ) c?ty^[P»(cost) + P»_i(cost)] sin refl. 

Since #)7r — e)t, equation ()2.14j) implies 

1 d f 9 r- £ hi(t)dt \ 

The evaluation of this integral is not immediate and is given in Appendix O 
We find that (ea. ljTHTjl ) 

/(a) = ^ l^/ (2^ + l)!) 2 J2^\ 

A } sVT^ (2^ + 2)! (2n+l)!^ } 

(2n)! (2n + 2)!(2n + 2) 2 \ ^n, nnW9 ^ 
" & jL ^(2^i)5 - ( 2 n + i(„ + i)!)2 a J ( X - a ) + ( 2 - 32 ) 

7T-6» . 

where a = , a < 1. The flux has a singular part, represented by 

e 

the half-integer powers of (1 — a 2 ), and a remaining regular part (the integer 

a 2 

powers.) The first term, , is the most singular one, because it 

eyl — a 2 

becomes infinite as \a\ — > 1. In other words, the flux is infinitely large near the 
boundary of the hole. The splitting of the solution into singular and regular 
parts is common in the theory of elliptic boundary value problems in domains 
with corners (see e.g., [2~0]-|22]). 

The value of the flux at the center of the hole is to leading order 



»-S(?H-^)^A ,,33) 



e ^ v. 2 (2 n n\) 2 (2n + l)!/ : 



in agreement with ()2.29|) (thanks Maple for calculating the infinite sum.) 
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The size of the boundary layer is varying with 9 proportionally to 1 //(#). 
The singularity at the end points of the hole indicate that the layer shrinks 
there to zero. Therefore, the boundary layer is shaped as a small cap bounded 
by the absorbing arc and (more or less) the curve y/1 — a 2 (see Fig|2J). In 

particular, the MFPT on the reflecting boundary is O ^log -j , even when 

taken arbitrarily close to the absorbing boundary. The singularity of the flux 
near the endpoints indicates that the diffusive particle prefers to exit near the 
endpoints rather than through the center of the hole. 

The expansion (jC.17j) is useful in approximating the flux near the endpoints 
(a = ±1), where few terms are needed. However, it is slowly converging near 
the center of the hole, where a power series in a 2 should be used instead 

oo 

/(«) = 5>a 2n + 0(l), (2.34) 

n=0 

where the coefficients f n are O^e" 1 ). Equations ()2.33|) and ()2.29|) indicate that 
/o = — . All other coefficients can be found in a similar fashion. We conclude 

£ 

that near the center (a <C 1) we have 

/(a) = ~ + o(l,y). (2.35) 

A Maximal exit time for the circular disk 

Using equation (|2.1()j) we find 



1 n =l 

= ^- + A= I hx{t) YV n (cost) +P n _ 1 (cost)]dt. 
1 v 2 Jo n=1 

Recall the generating function of the Legendre polynomials [T9] 

1 oo 

Y,Pn{x)t n , (A.2) 



oo 



VI - 2tx + t 2 

v n=0 



from which it follows that 



^P n (cost) = = p (A.3) 

^ Vl-2cost + l 2in l 

2 
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Together with equation (|2.11|) . this gives 



v max = — + -^=J^ hi{t) 



\ Sm 2 



\ 



/ 



dt = —= . 

V2J0 



hi(t) dt 



sin • 



(A.4) 



Combining with equation ()2.16j) and integrating by parts, we get 



\[2 Jo 71 a ;„ 1 dt J y/cos u - cos t 
8 2 

t u sin — du t-£ 
2 



dt 



tc sin ■ 



£ J \J cos u — COS t 



(A.5) 



u 



2V2 



7T Jo 



TVS COS ■ 



sin 



u sin — du 



jdt . 

_ Jo V cosu ~ cos t 



Equations (l2~TTI) and (l2~T31) show that 



u 



7T Jq y/cOSU — COSt 



2 log cos — h 2 log 1 + sin- )+k{t), (A.6) 



where 



k(t) 



sin 



/ 



Therefore, 



arcsm -\ / s 2 + cos 2 | 



\ 



s 2 + cos 2 - 



t 

2 / 



ds. 



(A.7) 



t u sin — <iu 
2 

7T sin - ^° Vcosm - cost 



lim 



v/2 



-2 log cos - + 2 log 1 + sin - + 



lim 

t->o 



sm 



2 arcsin(l) = 0. 

7T 



(A.f 
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Hence 



£ 

2 cos- 



2 log C 1 + cos ^) — 2 log sin ^ — — / 
\ 2/ 2 7r y 



cos- 



fvrcsin , is 2 + sin 2 



s 2 + sin 2 - 



t . u 

1 r-K-e COS — rt U Sin — GSM 



2\/27r Jo gin 2 ^ 7 Vcos - " — cost 
2 



2 / 



For e < 1 



-log 



1 



t . u . 

7T cos - r t u sin — aw 

2 - dt ^ 2 



2 2 v / 2r; ./,, gin 2 _ ./ 
2 



* 7o V cos u — COS t 



+ 0(e). (A.9) 



Changing the order of integration, we get 



, e i r 

v - = - log 2 + ^vrJo 



t , 

u r cos 2 dt 

u sin — du t 

2 / . o t / 

Ju sin -ycos-u — cost 



+ 0(e). (A.10) 



Substituting 



in the inner integral results in 



cosm — cost 



(A.ll) 



/ 

J u 



7T cos - dt 

2 



u 

cos — 
= V2- 2 



u ■ 



sin 2 -a/cosm — cost si n 77 



Therefore, 



i r 

^ Jo 



d u = _log_ - _ / 

2 27T ./ n +QT1 _ 2 7T J Q 



tan — 
2 



log - + log 2. 



log sin v dv 

(A.12) 
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B Exit times along the ray 

Along the ray 9 = tt the MFPT is given by 

1 r 2 °° 
v ray (r) = v(r,9 = n) = + ^ + ^2 a n(~r) n 

n=l 

-1 2 "1 pTT — £ 

= + v + ^iWE[ p «( cost )+ p -i( cost )K-o n ^ 

Using the generating function (|A.2|I of the Legendre polynomials to sum the 
infinite series, we obtain 

Vray{r) = ~~ + —=- / 75== , ■ 2 = " R1 

4 v2 Jo vl + 2r cost + r 2 

Combining with equation (|2.1fij) . integrating by parts, and hanging the order 
of integration gives 

1 — r 2 1 — r 

v ra y(r) = —— + — j ===== a 

4 2vl — 2r cose + r 2 

rfl-r) r~ £ . u , r _£ sin tdt 

it sin — aw 



\/27r 7o 2 ,/ u (1 + 2r cost + r 2 ) 3 / 2 >/costi — cost 



The substitutions s = yj cos u — cos t and x = v 2r s lead to 



sin t (it 2y / o)s~u+~cose 



1 + 2r cost + r 2 ) 3 / 2 v / cosu — cost (1 + 2r cosm + r 2 )Vl — 2r cose + r 2 
which implies that 

v ray(r) = ] -^- + —j= = a (B.2) 

4 2a/1 — 2r cose + r 2 

\/2r(l — r) Z^ -5 . u J cos u + cos e , 

tt sin att. 



7r\/l - 2r cose + r 2 Jo 2 1 + 2r cos u + r 2 

The substitution (j2.18j) gives 



tt a/cos tt + cos e 

ti sm att 

2 1 + 2r cos tt + H 



4^ 



cos £ arccos y s 2 + sin 2 ^ s 2 (is 



" (1 — 2r cose + r 2 + 4rs 2 ) ^ /sin 2 - + s 2 
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and we obtain the exact form of v ray (r) as 

v ra y{r) = ^—^ + = ===== a (B.3) 

4 2v 1 — 2r cose + r 2 



( , )S £ arccos ^s 2 + sin 2 - s 2 rfs 

7r\/l - 2r cose + r 2 J /, 2 , , 2 \ / • 2 5 , 2 



r (l_ r ) ^2 V 2 



(1 — 2r cose + r 2 + 4rs 2 )^/sin — + 



For e < 1 and 1 — r ^> yfe equation ()B.3|) becomes 

. . 1 — r 2 , £ 8r f 1 s arccos s ds . , s 

WO = — - log 5 - - jf (1 _. r)2 + 4rs , +0(e). (B.4) 



To evaluate the integral in ()B.4|) . we write 



7T 

arccos s = arcsin s, (B.5) 

2 



and obtain 



8r n f 1 sds . , , . 9 . 

V 2 X (l-r)' + 4r.« = ~ 2 ' 0g(1 - r) + l0S(1 + " 

The function g(r), defined by 

8r f 1 arcsin s sds 
q{r) = VJ (l-r) 2 + 4r, 2 (R6) 

in the interval < r < 1, has the endpoint values 

q(0)=0, g(l)=lo g 2. (B.7) 

Therefore, 

Vray(r) = - log | + 2 log(l - r) + ] —^- - log(l + r 2 ) + g(r) + 0(e), (B.8) 
is the MFPT for e C 1 and 1 — r ^> y^i. In particular, 

Vcenter = V ray (0) = - log | + - + 0{s) , 

as asserted in (jl.2|) . 
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C Flux profile 



In this appendix we calculate the flux profile given by equation (|2.3ip . Sub- 
stituting equation (j2.16|) for h\ in equation ()2.31|) gives 



f(0) = - 



1 d 

2 + d9 



1 d 



cos ■ 



t u sin — du 
2 



vr dt J cos u — cos t 
V cos t — cos 9 



dt 



Integration by parts and changing the order of integration, we find that 

9 



cos- 



u sin — du 
2 



a/cos(7t — e) — cos 9 Jo \J cos u 4~ cos e 

sin t dt 



1 
— cos — 

2 2./n -~"-~2""~J u (cost -cos#) 3 / 2 (cosu- cost) 1 ^ 



it sin — rfw 



We evaluate the inner integral by making the substitution x = ^cosu — cost, 
~ e smtdt 2-v/cosw + cose 



(cost - cos 9) 3 / 2 (cos it - cost) 1 ^ (cosu - cos0)V- cos6> - cose 



Therefore 



1 1 d 

2 + 

cos 



cos- 



m sin — g?w 
2 



V— cos £ — COS 9 Jo ^COS U + COS £ 



8 



u 



V cos — cose Jo 



7T-e u sin — a/ cos u + cos e 



cos m — cos 9 



du 



1 1 d 

2 + 

cos- 



9 

cos- 



7T 



V— cose — cos 9 \[2 



u 



a - 



V— cos 9 — cose 



tt-e u sin — i/cosm + cose 



cos it — cos 9 



du 
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The substitution (|2.18|) gives 



u 



TT~e u sin — V cos u + cos e du 
cos u — cos 6 



2V2 



cos ■ 



arccos \l s 2 + sin 2 - s 2 ds 



/ 2 1 • 2 e 2 

+ sin — — cos 



s 2 + sin 2 



Therefore, the flux takes the form 



/(*) 



1 1 d 

~z + 



2 V27cd6 



e 

cos - 
2 



V— COS6 1 — cose 



/ 



cos(e/2) 



7ra 



which is rewritten as 



4 arccos * /s 2 + sin 2 - s 2 ds 
1 2 



s z + sin cos - 



sin - + s 



7 



/(*) 



1 1 d 

2 + 2^d9 



cos- 



7TOo — 4 



arccos \/s 2 + a 2 s 2 ds 



Va 2 + s 2 (s 2 + b 2 ) 



(C.l) 



where a = sin - and 2b 2 = — cos 9 — cos e. Writine 

2 6 



<f){a,s) 



arccos 



Vs^T 



2n 



n=0 



we find the Taylor coefficients 

, , , arccos a . . / arccos a 1 

0o a = , 02 a = - 3 + ~ „ R g 

a V 2a 2a 2 VI - a 2 

and so on. For all n > we find the asymptotic behavior 

>2„ ( " i ~ -^TT + O I — ) ilS (/ — • 0. 



a 



2n+l 



2n. 
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To see this, consider the Taylor expansions 



2n+l 



S\* \ ^ s 2m 

a) I ^ n a 2m 

m=0 



/ • x 2n+l 



arccosjaWl +(-)") = £ + £ « n a 2n+1 ( Jl + ' 



2. — /t— II/ 1 \ 

n=0 

- + Va a 2n+1 Vc m — 

n=0 m=0 



00 9™ 00 
7T S ZV 



= 2 + T,-^J2 c >" a2 

m=0 n=0 

where a n and c™ are (known) constants, and 



S \2 a^-l (2 n n!) 2 a 2n 



Therefore 

1 °° fQ-nM e 2n /-7T °° « 2m °° \ 

YK ' a ^ y ' (2 n n\) 2 a 2n \ 2 ^o 2m ^ " / v ; 

n=0 v ' \ m=0 n=0 / 

from which it follows that 

t M = ±( { - ir ^ + {a) )-^. (C.4) 

This shows that 

as asserted. The asymptotic behavior (jC.5|) of the coefficients 02n(a) can be 
used to estimate the integral in equation ()C.1|) . 

s 2 + b 2 = i>" (a) I + (c ' 6) 

n=0 J{ ' 

To extract the asymptotic behavior of the integral as b — > 0, we use the long 
division 

c 2n+2 n (_ i\n+ll.2n+2 

^(_l)i^ a ^ + < V b (C.7) 



s 2 + fe 2 > s 2 + b 2 

j=0 
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and integrate it to yield 



Vl-a 2 g 2n+2 



+ b 2 



E 

3=0 



(C. 



n+li2n+2 



-l) n+L b 



^ ds 



S 2 + b 2 



E 

3=0 L 



2n-2j+l 



' 2n - 2j + 1 



+ (-l) n+1 6 2n+1 arctan 



The Taylor expansion 



arctan 



vT 



r 7r 



VI v — -i 

9+E 



/_]\m+l ^2m+l 



\2m+l 



b 2 ^ 2m + 1 (vT^a2) ; 

m=0 \ v ' 

gives the Taylor expansion of the integral (|C.8J) in powers of b as 



(C.9) 



Vl-a 2 s 2n+2 



S 2 + b 2 



ds = ^{-l) j b 2j - 



3=0 



2n - 2 j + 1 



+(-l) n+i 6 



n+li2n+l 



7T + ^(-l) m+l /r"' + l 



m=0 



2m + 1 (VT^a 



2A2m+l 



f\/l 7i2\2n-2j+l 

Y(-iy {V1 ~ a) — 

U 2n-2j + l 

oo 

+ (_l)«+l!L & 2n + l + ^ 



:-i) 



ri+m 



r2m+2n+2 



m=0 



(2m + l)(\/T r a f ) 2m+1 



Therefore, the Taylor expansion of the integral (jC.6|) is 

oo 

^2<p2n(a 



VT ^Ma,s]s 2 ds 



s 2 + b 2 



n=0 



+(-i) n+1 ^ 2 " +1 + E 



E(-d j 

L 3=0 



2\2n-2j+l 



,„ j0 (2m + l)(v / T^) 2m+1 
Rearranging in powers of b, we find that 



2n -2] + 1 

L2m+2n+2 







fa, s) s 2 ds 



n=0 



(CIO) 
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where the first three coefficients are 



Po(a) = > 02n(a) tt— =/ q>\a, s) ds = - a , 

n=o zn+L J 4 



Pi (a) 



7T0 o ( a ) 7rarccosa 



2 2a 



A(o) = _ E , 2 .„ (a) (^!r! + **») 



n=l 



2n- 1 

(f)(a, s) - (f) (a) ^ + o (a) 



s 2 v 7 !^ 
and all other coefficients f3 n are recovered in a similar fashion, 

7T 7T 2 f2lV 

^ = (- 1 ) J+1 2^^ = -T ( 2 4^ +Q(a ' 2J) ' 



fly 



-l)ME^( fl ) 2n-2i + l 



1 



i-i 

-^02n(a) 
n=0 

( " 1)5 (/ 



n=0 

1 

n=0 (2i-2n-l)(v / T^)2^2n-i 
^0(«,s)-El= o 02„(a) S 2 " 



1 



(2j-2n-l)(v / T^) 2 ^-\ 

We see that extra effort should be put in finding the even coefficients /3 2n - 
Expanding 

= ! ' -yWC + (en) 

rv ' ; 2^ Z-/( 2 n n !)2 2n + l ' V ; 

v n=0 

and noting that the following infinite sum has a regular contribution 

lim J- V (2n)! ( s2 + fl2 ) W da - tf. (C 12) 
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where Cj are constants (also can be written in term of hypergeometric func- 
tions), we find an alternative representation for the even coefficients, 



[a 



S ) ~ En=0^2n(a)s 2 



n=0 



(2j - 2n - lXv 7 ! 3 ^) 2 -?- 2 ™- 1 



_iy _Q + 0(a) + 



1 ^ (2n)! ( g 2 + a 2 )" ^ 

v^2 v /i2 T -T Z^ (2 n n! )2 2n + l ^W°J 8 

v n=0 v ' n=0 



s 2j 



2n 

ds 



j ~ l 1 

£^ 2nl ; (2j-2n-f)(v / T^) 2 ^ 2 "- 1 



vr 1 £1 (2n)! ( S 2 + a 2 )« £i , , 2 „ 



, (2W 2n + l -E^) 8 



-vCT 2 Vs 2 + a 2 ^(2"n!) 2 2n + . 
')'( / — "=° s2j ds 



o 



J 2(2i- 1 (j-l)!) 2 a 2 i- 1 



2j-2 



The integrals are given in [23J, 

ds i i=i (-i )» a - 1\ / s 2 \ n ^ +1/2 



s 2 Vs 2 + a 2 a 2 i^2n-2j + lV n yV^ + a 2 



nlj Z-^d On _ 9<7 4- 1 \ n \ s 2 < n 2 ' ^ ^ 



The binomial expansion gives 



f il±^l± = y 1 ( f\ p-Wa^ (C.14) 

Altogether, we find that the integral term in equation (j(J.f 3j) is 



I 



1^ y (2n)! ( g 2 + a 2 r ^ 

v^2 v ^ T ^ ^(2"n!) 2 2n + l 

= ^ = rfs 



2a 2 ^ '^2n-2j + l V n / V^- 1 
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This sum has the closed form [23j 



i=0 



2i + 1 V i 



[2 k k\) 2 



{2k + 1)1 



(C.15) 



and we have obtained the asymptotic form of the even coefficients 



02 



vr 1 (2^(j - l)!) 2 „ / 1 



O 



(C.16) 



2j 2a?i (2j - 1)! ' " V^" 1 
We are now able to find the asymptotic expansion of the flux profile flC.l|) . 



1 1 d 
~2 + 2^d9 



1 2 d 

2 ~ nd9 







cos • 



nao — 4 



arccos \/ s 2 + a 2 s 2 ds 



Va 2 + s 2 (s 2 + b 2 ) 



n=0 



Setting ea = tt — 9, we obtain after some manipulations that to leading order 
in small e the flux is given in the interval — 1 < a < 1 by 



/(«) 



a 



r^T^ 2 £ 



OO 



n=0 



(2" +1 (to+1)!) ; 
(2n + 2)! 



9 (2 n n! 
■or - 



'1 - a 2 r +1/2 



(2n)! (2n + 2)!(2n + 2) 2 

-a 



2e^[{2 n n\) 2 (2 n + 1 {n + l)!) 5 



(2n + l)! 
l-a 2 ) n + 0(l). (C.17) 
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Figure 1: A circular disk of radius R. The arclength of the absorbing boundary 
(dashed line) is 2eR. The solid line indicates the reflecting boundary. 
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Figure 2: The boundary layer, indicated by "BL" , is the area bounded by the 
absorbing boundary (dashed line) and the solid arc. Outside the boundary 

layer the MFPT is O (log - J . 
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